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Resumen. Los metodos numericos: mimeticos, diferencias finitas y elementos finitos, son analizados 
desde un punto de vista numerico. Se busca concluir sobre la eficiencia, orden de convergencia y costo 
computacional de estos metodos. El analisis es hecho en problemas de valor de frontera unidimensio- 
nales (ecuacion conveccion-difusion en regimen estacionario) con variaciones distintas en el gradiente, 
coeficiente difusivo y velocidad convectiva. 
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Metodos conservatives, Convergencia. 

Abstract. Numerical methods: mimetic finite differences and finite elements, are analyzed from a 
numerical point of view. It seeks to conclude on the efficiency, order of convergence and computational 
cost of these methods. The analysis is done in boundary value problems one-dimensional (convection- 
diffusion equation at steady) with different variations in the gradient, diffusion coefficient and convective 
velocity. 

Keywords: Mimetics methods, Finite Element methods, Finite differences methods, Conservative 
methods, Convergence. 


1 INTRODUCTION 

Los metodos numericos comunmente usados para definir aproximaciones numericas a la 
solucion de los diversos problemas de la ingenierfa y las ciencias, que se modelan a partir 
de ecuaciones diferenciales parciales son: el metodo de diferencias finitas (DF) y el meto¬ 
do de los elementos finitos (MEF). No obstante, en las ultimas dos decadas, un nuevo tipo 
de esquemas conservatives en diferencias finitas, conocido originalmente como operadores 
de soporte y posteriormente como metodos mimeticos (MIM) (Shashkov and Steinberg, 1995; 
Hyman and Shashkov, 1998; Hyman et al., 2002; Castillo and Grone, 2003), a mostrado su su- 
perioridad ante los esquemas clasicos de diferencias finitas (Freites, 2004; Guevara et al., 2005). 
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Cada uno de estos metodos presenta ventajas para algun tipo de problema y desventajas en 
otros. A pesar de esto, y sin introducir ningun tipo de sesgo, el MEF ha predominado en ca- 
si todas las aplicaciones que resultan de interes cientffico hoy en dfa. Y, en aquellas donde el 
metodo DF mantuvo su preponderancia, el MEF estandar ha sufrido modificaciones para op- 
timizar su aplicabilidad en estos problemas. Por ejemplo, se puede citar a (Fi, 2006; Riviere, 
2008) para problemas de dinamica de fluidos y (Cordero and Dfez, 2010) en el caso de proble¬ 
mas con conductividad discontinua. Por otro lado, mucho se ha dicho sobre la eficiencia de los 
metodos mimeticos en comparacion con el metodo DF (Arteaga and Guevara, 2008; Guevara, 
2005), principalmente en problemas difusivos estacionarios. Sin embargo, se ha dejado de lado 
su comparacion con el MEF, justificando el hecho a una alta complejidad de la implementacion 
computacional del mismo. 

En este trabajo, buscando una mejor interpretacion de las diferencias que existen entre los 
tres metodos, se presenta un analisis numerico de las propiedades que posee cada uno de estos 
metodos: convergencia, precision en la frontera del dominio, eficiencia en diversos problemas, 
flexibilidad al dominio ffsico, implementacion y costo computacional. Dicho analisis se enmar- 
ca dentro de la ecuacion de conveccion-difusion estacionaria para problemas unidimensionales; 
se busca facilitar el analisis de los metodos, dejando de lado, momentaneamente, las dificultades 
analfticas y geometricas que implica un dominio ffsico 2D. 

El resto del artfculo se estructura de la siguiente manera. En el siguiente apartado se intro¬ 
duce el problema de valor de frontera usado junto a sus condiciones de frontera (Robin). En el 
tercer apartado se describen brevemente los metodos numericos; posteriormente, en la cuarta 
seccion, se realiaza la experimentacidn numerica y la discusion de los resultados. Finalmente, 
se presentan las conclusiones obtenidas y las referencias del trabajo. 

2 PROBLEMA MODELO 

La ecuacion que modela la transferencia de calor por conveccion-difusion en regimen esta- 
cionario en su forma mas simple (problema unidimensional) viene dada por 

d 2 d 

k(x)—zu + u(x)—u = /, en £l = (a,b), (1) 

aar dx 

donde u(x) representa la temperatura (variable del problema) en un punto x del dominio Q = 
(a, b); k > 0, el coeficiente de difusion termica; v, la velocidad convectiva (o advectiva) y /, 
una funcion escalar que describe la existencia de una fuente o sumidero en el problema. La 
ecuacion se completa al definir las condiciones de frontera: 

a a u(a) + f3 a -^-u(a) = 7 a a b u(b ) + Pb-^u(b) = (2) 

con (ccj, /3i, 7 i}, i = a, b, parametros reales conocidos, y dependiendo de su valor se tendran 
condiciones de contorno Dirichlet, Neumann o Robin. Aunque la solucion general del problema 
de valor de frontera (l)-(2) se puede obtener en gran parte de las configuraciones del problema 
mediante metodos analfticos, el mismo resultara apropiado para dejar ver las propiedades de los 
metodos numericos en estudio, y simplificar la visualizacion y analisis de los resultados. 
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3 METODOS NUMERICOS 


En este apartado se describen brevemente los tres metodos utilizados. Una descripcion de- 
tallada del esquema mimetico a ser utilizado en este artfculo puede ser vista, entre otras, en 
(Guevara, 2005). Para el caso del MEF existe un sin numero de referencias que introducen el 
metodo desde distintas vertientes; aqui, nos inclinamos por (Becker et al., 1981; Solin, 2006; 
Calderon and Gallo, 2011). Por ultimo, en el caso de DF se remite al lector a (Strikwerda, 2004). 

3.1 Metodo Mimetico 

Los metodos mimeticos se basan en la discretizacion de los operadores clasicos de las EDP 
(divergencia, gradiente y rotacional) de tal forma que ellos satisfagan una version discreta del 
Teorema de Stokes o identidad de Green: 

(Dv,f) Q + {v,Gf) P = {Bv,f) I . (3) 

Aquf, D, G y B son las versiones discretas de sus continuos correspondientes: gradiente (V), 
divergencia (V-) y operador de frontera d/dn. Los () representan un producto interior genera- 
lizado con pesos Q, P y I. Usando la identidad (3) se obtiene una relacion para el operador de 
frontera 

B = QD + G'P. (4) 

Para la discretizacion, se define una malla cuya geometrfa esta dada por los nodos ay, con 
i = 0,1,..., N, y las celdas de la misma son los intervalos [ay_i, ay], para i = 1,..., N. 
El tamano de la celda, h, viene dado por h = 1/N, suponiendo que la malla esta distribuida 
uniformemente. Los nodos intermedios de las celdas queda dado por ay + i/ 2 = (ay+ay+i)/2, ver 
Figura 1. La solucion y el operador divergencia se definen en el centra de las celdas, mientras 
el operador gradiente en los nodos ay que definen las celdas (ver Figura 1). 
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Figura 1: Malla unidimensional mimetica. 


En este trabajo, se analiza el metodo mimetico que surge de los operadores discretos de se- 
gundo orden (tanto en frontera como en los nodos intemos) introducidos por (Castillo and Grone, 
2003) y dados en (5) y (6). 
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En (6) aparecen dos filas de ceros (una en la parte superior y otra en la parte inferior), con 
esto se busca obtener una matriz cuadrada al componer operadores. El operador de frontera B 
queda dado explfcitamente por 


B = 
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N+2xN+2 

A partir de las discretizaciones de los operadores, la aproximacion mimetica para la ecuacion 
de conveccion-difusion (1) queda dada por 

kBGU + vGU = (fcDG + vG)U = F, (7) 


donde U representa la solucion aproximada mimetica a la solucion exacta, u, del problema 

U = (U(x 0 ), U(x 1/2 ),..., U(x N _ 1/2 ), U{x N ))\ 

y F representa la restriccion de / a la malla mimetica: 

F = (f(x 0 ), f(x i /2 ), • • •, f{x N - i /2 ), f(x N ))\ 


Como el operador divergencia discretizado no actiia sobre la frontera, entonces las condicio- 
nes de contorno Robin (2) se obtienen a partir 

[[a]+ [/?](BG)] £/ = /„, (8) 

donde el vector f b resulta de restringir el termino no homogeneo de las condiciones de contomo 
a la malla mimetica, es decir, f b = (7 a ,0,... ,0,7;,)*. Las matrices [a] y \/3\ son tales que 
a l)X = a a , a N+2 ,N+2 = ce b , /3 1A = (3 a , /3 n+2 ,n+2 = A,, y el resto de entradas son cero. 

A partir de (7) y (8), el esquema mimetico para la ecuacion de conveccion-difusion (1) sujeto 
a las condiciones de contorno Robin (2) queda dado por 


a 


+ [/3](BG) + kBG + vG U = F + f b . 


(9) 
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3.2 Metodo de Diferencias Finitas 


En el metodo de diferencias finitas los operadores diferenciales presentes en la ecuacion 
diferencial son discretizados (aproximados) para los nodos de la malla, Xi, con i — 0,1,..., N, 
usando esquemas en diferencia que provienen del truncamiento del polinomio de Taylor. Para 
obtener el mismo orden de convergencia que el esquema mimetico (7) planteado en el apartado 
anterior, se propone un metodo usando los esquemas de diferencia de segundo orden: 

jf_ u(xi) = + d_ u(xi) = -«(*,-.) +0(ft ^ 

En los nodos interiores de la malla, la discretizacion realizada a partir del esquema de segundo 
orden es riguroso. No obstante, al usar el esquema de diferencias centradas de segundo orden 
para discretizar las condiciones de frontera de Robin (2) se emplean nodos fantasmas, lo cual 
hace al metodo poco riguroso en este respecto. 

El esquema en diferencias finitas para la ecuacion de conveccion-difusion (1) queda dado 
entonces 


h 2hot a 

(■«o - 5-2fc 0 

^ Pa 


Uo + 2koUi — h 2 fo + (k® — VQ—)—p a , 


( k i ~ Vi—)Ui-i — 2 kiUi + (hi + Ui~)Ui + 1 


2 /3 a 

= h 2 fi, para * = 1 : N — 1, (10) 
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donde Ui, con i — 0,1,..., N, define la solucion aproximada del problema. Los subfndices en 
el esquema denotan la evaluacion en los nodos x, de la malla de diferencias finitas (por ejemplo, 

fi = 


3.3 Metodo de los Elementos Finitos 

En el MEF se tiene una solucion u H := con a* coeficientes a determinar, y fa 

funciones base de un espacio de dimension finita % C Y, con Y el espacio de funciones 
admisibles de la forma debil o variacional del problema 

B(u,v) = l(v), Vw G Y, (11) 


donde 

B(u , v) 

es una forma bilineal y 


f b 

9 a ot a u{a ) — 6 b a b u(b ) — / [( kv)' + uv] u'dx 

J a 


K v ) = / f vdx - ®bib + ObPb 

J a 


( 12 ) 

(13) 


es un funcional lineal. En (12) y (13), 9 a y 9 b quedan dadas por 9 a = k{a)v{a)/(3 a y 9 b = 
k(b)v(b)//3 b , y Y el espacio de Sobolev H 1 (fl). Los coeficientes a, son determinados a partir 
del sistema de ecuaciones que surge al sustituir u H en la forma variacional (11). 
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4 RESULTADOS NUMERICOS 


En este apartado, se presentan y analizan los resultados numericos obtenidos (comparati- 
vamente) por los metodos propuestos en distintas configuraciones del problemas (l)-(2) en el 
dominio f) = (0,1). Para el analisis se usan las normas del maximo, || • ||oo> y la norma L 2 
definidas por 

2 do. 

En la norma del maximo, Uj representa la solucion obtenida por alguno de los tres metodos 
numericos y u 3 la solucion exacta del problema. Para el caso del esquema mimetico, j = i+ 1/2, 
con i — 0,... ,n — 1; y j = 1,..., n — 1, para los otros dos metodos numericos. Para la norma 
L 2 , se usa la forma continua en lugar de su version discreta, ||e|| 2 2 = h |Uj — Uj | 2 , usada por 

(Batista and Castillo, 2009). Se quiere evitar el sesgo que puede presentar la version discreta de 
la norma a favor del esquema mimetico o diferencias Anitas. Pues, en cada segmento de la malla, 
la solucion de estos metodos es constante, mientras la de elementos finitos, no necesariamente. 

Ejemplo 1 

Se resuelve (l)-(2) con coeficiente de difusion k — 1, velocidad convectiva u — 0, condicio- 
nes de contorno Robin: 

cra(0) - i/(0) = —20/(e 20 - 1), au{ 1) + m'(1) = 0, 

con a := —20e 20 /(e 20 — 1) y un termino fuente, f(x), definido tal que la solucion analftica 
del problema queda dada por u(x ) = (e 20x — 1)/(e 20 — 1). La Figura 2 muestra las soluciones 
aproximadas junto a la solucion analftica del problema en una malla uniforme de 20 puntos. La 




Figura 2: Ejemplo 1: Comparacion de las soluciones aproximadas en un malla de 20 nodos. En la frontera, el FEM 
supera en precision a los otros dos esquemas. 


ampliacion en la frontera derecha muestra la mejor precision del MEF al imponer condiciones 
de frontera, seguido del esquema mimetico que supera a DF, aunque en esta se este usando es¬ 
quemas de segundo orden para aproximar la frontera (ver (10)). Esta propiedad de aproximacion 
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se mantiene al aumentar el numcro de nodos y se repite para cualquier problema con coeficien- 
tes k y v constantes. La igualdad en la precision de los nodos de frontera es solo alcanzada 
cuando se tienen condiciones de frontera Dirichlet. 

La Figura 3 muestra los errores alcanzados (a la izquierda medidos a partir de la norma del 
maximo y, a la derecha, usando norma L 2 ). En la grafica, las pendientes de las rectas definen 
el orden de convergencia de los metodos y mientras mas abajo se encuentre, mejor sera el 
esquema que represente. Los metodos alcanzan el orden de convergencia teorico en la norma 




Figura 3: Ejemplo 1: Errores numericos: en noma del maximo, || • Ijoo, (izquierda) y norma L 2 (derecha). 


N° de 

Elementos 

Norma • | oo 

Norma || ■ ||l 2 

DF 

MIM 

MEF 

DF 

MIM 

MEF 

1000 

1.999944 

2.016795 

4.000306 

0.986027 

1.000097 

1.999880 

3000 

1.999990 

2.007299 

4.022705 

0.994014 

1.000020 

1.999978 

5000 

1.999998 

2.003188 

3.485599 

0.997404 

1.000004 

1.999996 

7000 

2.000000 

2.002136 

-10.558 

0.998307 

1.000002 

1.999967 

9000 

1.999999 

2.001412 

15.251 

0.998740 

1.000001 

2.000043 

11000 

1.999998 

2.001592 

-27.111 

0.998995 

1.000001 

2.000664 


Cuadro 1: Ejemplo 1. Orden de convergencia en Norma del Maximo y Norma L 2 . 


del maximo (ver Tabla 1). Sin embargo, en norma del maximo, para este problema, el MEF 
presenta una superconvergencia alcanzando el e de maquina. Por tal motivo, las oscilaciones 
que presenta la grafica para el caso del MEF es debido a errores de redondeo de maquina y 
no a perdida de precision del MEF. Esta caracterfstica no se mantiene en otras configuraciones 
del problema como se vera posteriormente. En norma L 2 , el MEF presenta un orden dos en su 
convergencia (como era de esperar, ver (Solfn, 2006)), mientras los otros dos metodos presentan 
un orden uno de convergencia (Figura 3 y Tabla 1). 

Ejemplo 2 

En este caso, se define el problema modelo (l)-(2) con coeficiente de difusion k(x) = 1/a + 
a(x — x 0 ) 2 , velocidad convectiva is(x) = k'{x), con a = 250 y x 0 = 0.75, condiciones de 
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contomo Robin: 


w(0) + u\ 0) = a/(l + a 2 Xo), w(l) + t/(l) = — arctan(a(l — tc 0 )) — arctan(a'X 0 ), 

y un termino fuente, f(x), definido tal que la solucion analftica del problema viene dada por 
u(x) = (1 — x) [arctan(a(x — x 0 )) + arctan(cta; 0 )] • La Figura 4 (izquierda) muestra las so- 
luciones aproximadas junto a la solucion analitica del problema para una malla uniforme de 60 
puntos. A la derecha, se ilustra la convergencia asintotica, h —>■ 0, de las aproximaciones en el 
nodo de frontera x — 1. Cuando la discretizacion es gruesa (pocos nodos) el metodo mimetico 
supera levemente al MEF. El metodo DF necesita superar los 1500 elementos para alcanzar la 
precision de los otros dos metodos. Para este ejemplo, se puede ver que, w(0) = w(l) = 0, lo 



Figura 4: Ejemplo 2: Soluciones aproximadas y exacta en un malla de 60 nodos (izquierda). Convergencia de los 
metodos en el nodo de frontera x = 1 a medida que se aumenta el numero de nodos en la malla. 

cual simplifica las condiciones de frontera y, por lo tanto, la formulacion variacional del pro¬ 
blema (ver, ecuaciones (12)-(13)). En estos casos, la mejor precision en los nodos de frontera 
es lograda por el MEF 

La Figura 5 y Tabla 2 muestra los errores y el orden de convergencia alcanzado. En norma del 




Figura 5: Ejemplo 2: Errores numericos: en noma del maximo, || • ||oo, (izquierda) y norma L 2 (derecha). 


maximo, el metodo DF pierde el orden 2 de convergencia de los esquemas usados para su cons- 
truccion. Este fenomeno es debido a que los coeficientes de difusion y convectivo dependen de 
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la variable espacial (no constantes). El MIM y MEF mantienen el orden 2 de convergencia con 
MEF, superando ligeramente en precision a MIM. En norma L 2 , se repite el comportamiento 
del ejemplo anterior: orden dos para MEF y orden uno para DF y MIM. 


N° de 

Elementos 

Norma • 

Norma || ■ \\l 2 

DF 

MIM 

MEF 

DF 

MIM 

MEF 

1000 

1.154080 

1.942273 

1.993329 

1.004321 

1.003528 

1.996938 

3000 

1.000175 

2.004354 

1.998233 

0.998419 

1.000636 

1.998873 

5000 

1.000073 

1.994971 

1.999533 

0.998652 

1.000109 

1.999803 

7000 

1.000047 

2.000917 

1.999993 

0.999021 

1.000045 

1.999919 

9000 

1.000035 

2.000485 

1.999783 

0.999236 

1.000025 

1.999950 

11000 

1.000028 

1.998615 

2.000046 

0.999375 

1.000016 

1.999978 


Cuadro 2: Ejemplo 2. Orden de convergencia en Norma del Maximo y Norma L 2 . 


Ejemplo 3 

Se resuelve el problema modelo (l)-(2) para un coeficiente de difusion constante k = 1.052, 
velocidad convectiva v = —110.5, y un termino fuente / = 0. La solucion analftica queda 
dada por u(x) = (1 — e~ Xx )/(l — e _A ), con A = vfk. En la experimentacion numerica se to¬ 
man inicialmente condiciones de contomo Dirichlet: ?j(0) = 0, u (1) = 1; y, posteriormente, 
condiciones Robin: 

u'( 0) = A/(l - e _A ), u( 1) + u\ 1) = 1 + \e~ x /(1 - e" A ) 

Para esta configuracidn del problema junto a condiciones de contomo Dirichlet, es altamente 
conocido y facil de probar, los metodos DF y MEF resultan equivalentes y presentan oscilacio- 
nes cuando el numero local de Peclet \v\h/(2k) 1. La forma mas simple de superar este 

fenomeno oscilatorio (ffsicamente no correcto) es hacer h suficientemente pequeno o anadien- 
do un termino de difusividad artificial al problema. Debe quedar claro, que en los casos en que 
el numero local de Peclet es igual o mucho menor que 1 estos metodos no presentan oscila- 
ciones, contrario a lo que erroneamente se afirma en (Batista and Castillo, 2009). La Figura 6 
muestra las soluciones obtenidas para 50, 80 y 200 elementos en la malla. Resulta evidente que 
la solucion del esquema mimetico tambien presenta las oscilaciones, y para este caso necesita 
mas nodos que DF y MEF para evitas las oscilaciones y lograr la convergencia a la solucion 
analftica. 

La Figura 7 y Tabla 3 muestran los errores y el orden de convergencia alcanzado para el 
caso de condiciones de contomo Dirichlet. En norma del maximo, como ya se dijo antes, los 
metodos DF y MEF resultan equivalentes, y mantienen el orden 2 de convergencia. Sin embargo, 
el esquema MIM pierde su orden 2 de convergencia y su precision es mas pobre que la lograda 
por los otros dos metodos. En norma L 2 , el metodo DF se une a MIM logrando solo un orden 
uno en su convergencia, y con un error exiguo entre los dos. En esta norma, el MEF mantiene 
el orden dos de convergencia y nuevamente resulta superior en su precision. 
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Figura 6: Ejemplo 3: Soluciones aproximadas y exacta para 50 (izquierda), 80 (centra) y 200 (derecha) elementos. 



Figura 7: Ejemplo 3: Errores numericos: en noma del maximo, || ||oo, (izquierda) y norma L 2 (derecha). 


N° de 

Elementos 

Norma || • ||oo 

Norma ■ \\l 2 

DF 

MIM 

MEF 

DF 

MIM 

MEF 

1000 

2.005097 

1.078854 

2.005097 

0.990024 

1.042730 

1.998561 

2000 

1.999493 

1.034747 

1.999493 

0.995134 

1.018963 

1.999640 

4000 

2.000324 

1.016534 

2.000324 

0.997599 

1.008977 

1.999910 

5000 

2.000267 

1.009977 

2.000267 

0.998517 

1.005461 

1.999966 


Cuadro 3: Ejemplo 3. Orden de convergencia en Norma del Maximo y Norma L 2 . 


Para las condiciones de contorno Robin, las soluciones aproximadas son mostradas en la 
Figura 8. En la frontera x = 1 la precision es pobre para los metodos DF y MIM (parte ampliada 
al centra de cada grafica de la Figura 8). El esquema MIM llega a necesita hasta 4 mil elementos 
para alcanzar la precision que logra MEF con solo 80 elementos. EL metodo DF pierde precision 
en esta frontera y no converge. 

Los resultados de convergencia para MEF y MIM resultan equivalentes a los alcanzados 
para el caso de condiciones de contorno Dirichlet (orden 2 para MEF y orden 1 para MIM). 
Sin embargo, diferencia Anitas pierde por completo su convergencia, debido a la perdida de 
precision que presenta alrededor de la frontera x — 1. Por ultimo, se debe senalar que la matriz 
que define el sistema resultante para MEF y MIM presenta un numero de condicion superior a 

10- 19 . 
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Figura 8: Ejemplo 3: Soluciones aproximadas y exacta para 100 (izquierda) y 200 (derecha) elementos. La parte 
ampliada en cada grafica representa el comportamiento numerico alrededor de la frontera x = 1. 

5 CONCLUSIONES Y COMENTARIOS FINALES 

Resulta claro, que el MEF presenta mejores resultados en precision y convergencia que los 
otros dos metodos. A1 mismo tiempo, MIM resulta en todo momento superior al metodo DF 
(como ya ha sido reportado por muchas referencias). En el caso del Ejemplo 3, donde se podrfa 
esperar superioridad del MIM, debido a su condicion conservativa, tampoco logra superar a 
MEF. Sin querer entrar en discusion de las modificaciones apropiadas o metodos optimos para 
la resolucion del Ejemplo 3, queda abierta la pregunta de que condiciones o cambios se deben 
imponer al MIM para problemas altamente convectivos. 

En gran medida, se ha justificado el uso de MIM ante MEF debido a la complejidad teorica 
de MEF (para aquellos no matematicos). Sin embargo, la definicion de los operadores discretos 
usados por MIM (principalmente en el caso 2D o 3D) puede resultar tan complejo en calculo 
que pueden competir facilmente con la necesidad teorica de MEF. Tal dificultad se acrecienta, 
si se quiere subir el orden de convergencia del metodo, hecho que se logra en MEF con gran 
facilidad. Ademas, hay que anadir la obligatoriedad que existe de una eficiente implementacion 
computacional, pues el uso de matrices para definir operadores y condiciones de contomo puede 
resultar en un rapido desbordamiento de memoria en el metodo MIM. Por ultimo, y tal vez 
lo mas importante, se debe mencionar el alto grado de complejidad y dificultad que presenta 
MIM para trabajar en mallas que se ajusten a dominios arbitrarios o en elementos deformados 
(Hyman and Shashkov, 1997, 1998). 

Como conclusion, si se quiere abordar un problema con un metodo novel y se tiene a mano 
la discretizacion de los operadores, el metodo mimetico, MIM, resultara una buena eleccion, 
por arriba de cualquier esquema de diferencias Anitas. Desde la vertiente matematica, sobresale 
el interes de afianzar las bases teoricas del metodo MIM y proponer mejoras en las lfneas donde 
el mismo aun sigue sin rigor o no han sido abordadas. 
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